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We report on the results obtained in the study of the parton saturation effects, 
taken into account through the multiscattering Glauber-Mueller approach, applied 
to the Drell-Yan (DY) process described in the color dipole picture. As a main 
result, one shows that those effects play an important role in the estimates of the 
DY differential cross section at RHIC energies. 



The Drell-Yan (DY) process, i.e. the production of massive lepton pairs 
in hadronic collisions, in conjunction with deep inelastic scattering (DIS), 
has been one of the important processes probing strong interaction physics. 
Recently, in connection with the availability of high energy accelerators, a 
great attention has been focused on the small- a; region of QCD. There, the 
parton densities become high and the limits of the perturbative methods are 
tested. Such a region presents the onset of the saturation phenomenon (at 
scale Qs), i.e. the taming of the parton (mostly gluon) distribution due to 
nonlinear dynamics associated with unitarity effectso . 

In this contribution, we study the high energy DY cross section in the tar- 
get rest frame. In this case, the relevant degrees of freedom are the projectile 
wavefunction and the dipolc-proton effective cross section. The underlying 
process is the scattering of a parton from the projectile structure function 
off the target color field. This parton radiates (bremsstrahlung) a massive 
photon, which subsequently decays into a lepton pair. The interaction with 
the target can occur before or after the photon emission. A remarkable fea- 
ture emerging is that the 7*g-proton (or hadron) interaction can be described 
by the same qq (color dipole)-proton cross section as in DIS. Although dia- 
gramatically no dipole is present, the interference among graphs results in a 
product of two quark amplitudes in the DY cross section testing the external 
gluonic field at twa different transverse positions (impact parameters), in a 
similar way to DIScl. 

In such a representation, the photoabsorbtion cross section on deep inelas- 
tic is described by the convolution of the wavefunctions, ^'^y* , from the virtual 
photon and the interaction dipole cross section, a dip- The wavefunctions are 
considered taking into account the simplest photon Fock state configuration, 
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i.e., a qq-paii fluctuation (dipole), whose transverse separation r± remains 
fixed during the interaction since its lifetime is larger than the collision one. 
The adip is modeled phcnomcnologically based on a matching between hard 
and soft domains, constrained by the DIS available data. Small dipole size 
configurations can be described through pQCD, whereas the large size ones 
belong to the nonperturbative domain. Hence, one can write the photoab- 
sorptions cross section as a function of the Bjorken scaling variable x — xsj, 
the photon vituality scale = and the quark momentum fraction a in 
the quantum mechanical formtl, 

<Jt,l{i*P-^ qq) = J d^r± J da \-^'^^^{a, r±, Q^)^ adip{x, r±) , (1) 

where the simplest phenamenological realization of the dipole cross section is 
given by the GBW modec, a^ip = crQ[l — exp{—r^Ql/4)], with the parameters 
determined by fitting small-x HERA data. 

In a similar way, the cross section for radiation of a virtual photon from a 
quark after scattering on Euproton, has the following factorized form in impact 
parameter representation eI, 

^ d rj_ |^'^.^(a,ri, Q )\ (Tdip{x, ar±), (2) 

where '^■yq is the quark wavefunction representing the j*q fluctuation. Here, 
r_L is the photon-quark transverse separation, ar± is the qq separation and a 
is the fraction of the light-cone momentum of the initial quark taken away by 
the photon. Note the difference with the DIS case, where the dipole separation 
is just r± and now the dipole cross section is evaluated for ar±. The scaling 
variable x is taken as X2 — xi — xp, which are the usual invariants in the 
fast proton system: Xi^2 are the longitudinal momentum fractions of the 
projectile (target) parton and xp is the Feynman variable. We stress that such 
identification is far from clear in the dipole picture: for instance, one can also 
choose X = ax2- In the rest frame, the quark carries momentum fraction xi/a 
of the parent hadron, where xi is the proton's momentum fraction carried 
by the foton. The hadronic diffecential cross section to Drell-Yan process is 
expressed in the factorized form B, 

da^^ _ ckom 1 f'^da^pfxi ^2\ d<y{qp ^ ^* X) 



dAP dxp 6t:M^ (xi + X2) J^^ 01 V a ' / dhia 

where the factor aom/(6 7rM^) is due to the photon decay into a lepton pair, 
coming from electrodynamics. The dilepton invariant mass is Af^. Our input 
for (Jdip is constrained by the standard gluon distribution from the target 
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corrected by saturation effects in high energy Umit, encoded in the Glauber- 
Mueller (GM) approach □. In Eq. (||), the structure of the projectile is 
described by the proton inclusive structure function F2{x,Q^), in which was 
chosen = . Other possible identification relying on plausive arguments 
is = (1 — xi)M^ 13. Our input for the dipole cross section is given by, 

a^^{i,r^)^^r'xG%'\£,^), (4) 

7^ Jx X Ji/Q^ rl 

where 7e and Ei (x) are the Euler constant and the exponential integral. The 
packing factor is — {3TTasr^/2R^) xG^'^^^^ . The saturation scale is 
defined through the relation kg(x,Q^) — 1. The target transverse size is R, 
which was taken toJpe — 5 GeV~^ (strong unitarity corrections) supported 
by previous studiesE . The xG was parametrized by the GRV94 gluon density. 

In a recent workB, we have calculated the DY cross section for the available 
small-X2 data, and found that the dipole picture in that kinematical region ( 
X2 > 10^^) is being tested in its limit of applicability. As an improvement, 
we have introduced a reggeonic piece to account for the large X2 range in the 
phenomcnological form (7jj(a;,rj_) = Nji 7-2^2;0-4525^-[^ _ ^.^s^ which supplement 
Eq.(H). Towards smaller X2, such contribution turns out to be negligible. In 
Ref.L[ the procedure was to correct the saturation scale by a threshold factor 

—> (1 — x)^ , and the numeric equivalence between the. dipole picture 
(using GBW model) and NLO DGLAP calculations is shownQ, relying on an 
educated guess for the scales x and Q^. ^ 

In order to address the color dipole picture at high energies^ (small X2), 
the DY differential cross section for RHIC energies, y/s — 500 GeV, is shown 
in Fig. (|l|) for the same fixed xp- In this case, X2 reaches values of order 10~^ 
and unitarity effects are important. The solid lines are the GM estimates, Eq. 
(^, and the dot-dashed lines are the rest frame calculations with DGLAP 
gluon distribution (without parton saturation effects). The GM results un- 
derestimates those ones without parton saturation, showing important role of 
the unitarity effects. Moreover, we have obtained a result numerically equiv- 
alent to the LO Breit frame calculations at the RHIC energies, suggesting 
consistency between both approaches. From the plots one verifies that the 
deviations due to saturation effects are larger as AI diminishes, correspond- 
ing to smaller X2 values. In absolute values, the corrections at RHIC energies 
reach up to ~ 20 %, showing that saturation effects are important in a reliable 
description of the observables. 
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Figure 1. The DY differential cross section per nucleon versus M for the RHIC energies 
{^/s = 500 GeV) at fixed xp in pd reaction. The solid line are GM estimates, whereas 
dot-dashed ones are without gluon saturation calculation. 
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